Particle showers detected on ISS in Timepix pixel detectors.

We detect regular particle showers in several compact pixel detectors, distributed over the International Space Station. These showers are caused by high energy galactic cosmic rays, with energies often in the 10 s of TeV or higher. We survey the frequency of these events, their dependence on location on ISS, and their independence of the location of ISS, on its orbit. The Timepix detectors used allow individual particle tracks to be resolved, providing a possibility to perform physical analysis of shower events, which we demonstrate. In terms of radiation dosimetry, these showers indicate certain possible limitations of traditional dosimetric measures, in that (a) the dose measured in small sensor may be less than that received in a larger distribution of matter, such as a human and (b) the spatial and temporal extent of these events represents a regime of poorly documented biological response.

Hybrid methods of radiation shielding against deep-space radiation.

In the last decade, NASA and other space exploration organizations have focused on making crewed missions to different locations in our solar system a priority. To ensure the crew members' safety in a harsh radiation environment outside the protection of the geomagnetic field and atmosphere, a robust radiation protection system needs to be in place. Passive shielding methods, which use mass shielding, are insufficient as a standalone means of radiation protection for long-term deep-space missions. Active shielding methods, which use electromagnetic fields to deflect charged particles, have the potential to be a solution that can be used along with passive shielding to make deep-space travel safer and more feasible. Past active shielding studies have demonstrated that substantial technological advances are required for active shielding to be a reality. However, active shielding has shown potential for near-future implementation when used to protect against solar energetic particles, which are less penetrating than galactic cosmic rays (GCRs). This study uses a novel approach to investigate the impacts of passive and active shielding for protection against extreme solar particle events (SPEs) and free-space GCR spectra under solar minimum and solar maximum conditions. Hybrid shielding configuration performance is assessed in terms of effective dose and radiobiological effectiveness (RBE)-weighted dose reduction. A novel electrostatic shielding configuration consisting of multiple charged planes and charged rods was chosen as the base active shielding configuration. After a rigorous optimization process, two hybrid shielding configurations were chosen based on their ability to reduce RBE-weighted dose and effective dose. For protection against the extreme SPE, a hybrid active-passive shielding configuration was chosen, where active shielding was placed outside of passive shielding. In the case of GCRs, to gain additional reduction compared to passive shielding, the passive shielding configuration was placed before the active shielding to intentionally fragment HZE ions to improve shielding performance.

Modeling a 3-D multiscale blood-flow and heat-transfer framework for realistic vascular systems.

Modeling of biological domains and simulation of biophysical processes occurring in them can help inform medical procedures. However, when considering complex domains such as large regions of the human body, the complexities of blood vessel branching and variation of blood vessel dimensions present a major modeling challenge. Here, we present a Voxelized Multi-Physics Simulation (VoM-PhyS) framework to simulate coupled heat transfer and fluid flow using a multi-scale voxel mesh on a biological domain obtained. In this framework, flow in larger blood vessels is modeled using the Hagen-Poiseuille equation for a one-dimensional flow coupled with a three-dimensional two-compartment porous media model for capillary circulation in tissue. The Dirac distribution function is used as Sphere of Influence (SoI) parameter to couple the one-dimensional and three-dimensional flow. This blood flow system is coupled with a heat transfer solver to provide a complete thermo-physiological simulation. The framework is demonstrated on a frog tongue and further analysis is conducted to study the effect of convective heat exchange between blood vessels and tissue, and the effect of SoI on simulation results.

Methods for estimating X-ray machine output through measurement and simulation.

Ball grid arrays are increasingly being applied in the electronics industry and may require X-ray inspection to ensure the integrity and correct placement of solder pins. However, as the architecture of integrated circuits continues to narrow while simultaneously growing more complex, the risk of electronic failure due to radiation damage increases. While medical X-ray devices have been held to high standards and are repeatedly shown to be well characterized, devices used for electronic inspection are often lacking detailed characterization. This study presents unique methods to solve for important properties in X-ray inspection devices such as source to object distance and energy spectrum. This information can then be applied to Monte Carlo models to achieve better overall dose estimates to electronics, which will lead to superior manufactured products. Since X-ray devices can vary greatly in source characteristics, this work investigates spectral measurement and Monte Carlo representation of three X-ray devices. For a Philips SRO 33 100 medical diagnostic device, the spectral output followed expected trends given by the prediction software SpekCalc and Spektr. For the Dage XD7500NT, direct measurement showed a spectral artifact that through the use of Gafchromic films, was shown to be a contributing effect in the dose output. For the Rad Source RS1800, a high powered irradiation device, direct spectral measurement was not achieved. However, a Monte Carlo model using an assumed spectra was found to match ion chamber measurements to a high degree.

Bulk material interrogation experimental results and validation with Geant4 for replacement of dangerous radiological sources in oil-well logging industries.

The Kansas State University Materials Interrogation (KSUMI) test facility was set up to enable bulk-material irradiation experiments that replicate similar oil-well logging scenarios, with an aim to address the problem of replacement of conventional radioisotope sources commonly used in oil-well logging industries. An exploration tool similar to an oil-well logging tool was used to conduct experiments with water and sand as testing materials. The facility includes a 2500-gallon concrete test chamber with an aluminum pipe going horizontally through it. A machine-based 14.1 MeV deuterium-tritium neutron source as an alternative to conventional neutron sources was used. High energy neutrons assist in the investigation of a larger volume of material and also generate high energy inelastic scatter gamma rays, which provide useful information on composition. Experiments were performed with tap water and sand as a bulk testing material. Irradiation was done for one hour and results were obtained from a 3He neutron sensor, a BF3 neutron sensor, and two NaI gamma sensors placed at different locations within the exploration tool. Geant4, a Monte-Carlo based toolkit, was deployed on a high-performance computing system to simulate the entire experiment in order to benchmark the experimental responses obtained from the photon and neutron sensors. The facility was modeled in detail with accurate dimensions and material compositions. Materials such as tap water, high-density polyethylene, and aluminum metal were modeled with thermal neutron scattering cross-sections. The reference physics list QGSP_BIC_HP along with G4NDL and S(α,ß) cross-sections were found to be appropriate for simulation of neutron interrogation experiments with neutron energies lower than 20 MeV. The experimental results obtained were successful in characterizing the bulk testing materials, and results obtained from Geant4 were found to be in good agreement with the experimental results in most cases.

A Structured Cleaving Mesh for Bioheat Transfer Application.

Goal: The thermoregulation mechanism is a complex system that executes vital processes in the human body. Various models have been proposed to simulate the thermoregulatory response of an adult human to environmental stimuli. However, these models generally rely on stylized phantoms that lack the anatomical details of voxel phantoms used in radiation dosimetry and shielding research. The goal of this work is to introduce voxel phantoms to thermoregulation research by modeling the physical energy exchange between tissue and its surroundings, discuss a specific challenge associated with voxel phantoms, propose a method to address this challenge, and demonstrate its application. Method: One of the major challenges in using voxel phantoms is the stair-step effect on the surface of the voxelized domain. This effect causes over-estimation of surface area, accurate knowledge of which is critical for modeling heat exchanging systems. A methodology to generate a voxel domain from medical imaging data and reduce error in the surface area caused by the stair-step effect is presented. The methodology, based on a structured mesh and finite-volume method, is demonstrated with tumors generated from magnetic resonance imaging (MRI) scans of mice. Results: The methodology discussed in the paper shows a decrease in surface area over-estimation from 50% to 15% for a sphere and 47% to 17% for tumor models generated directly from MRI scans. Conclusion: This work provides a direct method to generate a smoother domain from medical imaging data and reducing surface area error in a voxelized domain. The technique presented is independent of domain material, including tissue type, and can be extended to any homogeneous or inhomogeneous domain. The increase in surface area accuracy obtained by smoothing the voxel domain results in more accurate temperature estimates in heat transfer simulation.

Comparison of methods for individualized astronaut organ dosimetry: Morphometry-based phantom library versus body contour autoscaling of a reference phantom.

One of the hazards faced by space crew members in low-Earth orbit or in deep space is exposure to ionizing radiation. It has been shown previously that while differences in organ-specific and whole-body risk estimates due to body size variations are small for highly-penetrating galactic cosmic rays, large differences in these quantities can result from exposure to shorter-range trapped proton or solar particle event radiations. For this reason, it is desirable to use morphometrically accurate computational phantoms representing each astronaut for a risk analysis, especially in the case of a solar particle event. An algorithm was developed to automatically sculpt and scale the UF adult male and adult female hybrid reference phantom to the individual outer body contour of a given astronaut. This process begins with the creation of a laser-measured polygon mesh model of the astronaut's body contour. Using the auto-scaling program and selecting several anatomical landmarks, the UF adult male or female phantom is adjusted to match the laser-measured outer body contour of the astronaut. A dosimetry comparison study was conducted to compare the organ dose accuracy of both the autoscaled phantom and that based upon a height-weight matched phantom from the UF/NCI Computational Phantom Library. Monte Carlo methods were used to simulate the environment of the August 1972 and February 1956 solar particle events. Using a series of individual-specific voxel phantoms as a local benchmark standard, autoscaled phantom organ dose estimates were shown to provide a 1% and 10% improvement in organ dose accuracy for a population of females and males, respectively, as compared to organ doses derived from height-weight matched phantoms from the UF/NCI Computational Phantom Library. In addition, this slight improvement in organ dose accuracy from the autoscaled phantoms is accompanied by reduced computer storage requirements and a more rapid method for individualized phantom generation when compared to the UF/NCI Computational Phantom Library.

Optimal shielding thickness for galactic cosmic ray environments.

Models have been extensively used in the past to evaluate and develop material optimization and shield design strategies for astronauts exposed to galactic cosmic rays (GCR) on long duration missions. A persistent conclusion from many of these studies was that passive shielding strategies are inefficient at reducing astronaut exposure levels and the mass required to significantly reduce the exposure is infeasible, given launch and associated cost constraints. An important assumption of this paradigm is that adding shielding mass does not substantially increase astronaut exposure levels. Recent studies with HZETRN have suggested, however, that dose equivalent values actually increase beyond ∼20g/cm2 of aluminum shielding, primarily as a result of neutron build-up in the shielding geometry. In this work, various Monte Carlo (MC) codes and 3DHZETRN are evaluated in slab geometry to verify the existence of a local minimum in the dose equivalent versus aluminum thickness curve near 20g/cm2. The same codes are also evaluated in polyethylene shielding, where no local minimum is observed, to provide a comparison between the two materials. Results are presented so that the physical interactions driving build-up in dose equivalent values can be easily observed and explained. Variation of transport model results for light ions (Z ≤ 2) and neutron-induced target fragments, which contribute significantly to dose equivalent for thick shielding, is also highlighted and indicates that significant uncertainties are still present in the models for some particles. The 3DHZETRN code is then further evaluated over a range of related slab geometries to draw closer connection to more realistic scenarios. Future work will examine these related geometries in more detail.

Solar proton exposure of an ICRU sphere within a complex structure part II: Ray-trace geometry.

A computationally efficient 3DHZETRN code with enhanced neutron and light ion (Z ≤ 2) propagation was recently developed for complex, inhomogeneous shield geometry described by combinatorial objects. Comparisons were made between 3DHZETRN results and Monte Carlo (MC) simulations at locations within the combinatorial geometry, and it was shown that 3DHZETRN agrees with the MC codes to the extent they agree with each other. In the present report, the 3DHZETRN code is extended to enable analysis in ray-trace geometry. This latest extension enables the code to be used within current engineering design practices utilizing fully detailed vehicle and habitat geometries. Through convergence testing, it is shown that fidelity in an actual shield geometry can be maintained in the discrete ray-trace description by systematically increasing the number of discrete rays used. It is also shown that this fidelity is carried into transport procedures and resulting exposure quantities without sacrificing computational efficiency.

Solar proton exposure of an ICRU sphere within a complex structure Part I: Combinatorial geometry.

The 3DHZETRN code, with improved neutron and light ion (Z≤2) transport procedures, was recently developed and compared to Monte Carlo (MC) simulations using simplified spherical geometries. It was shown that 3DHZETRN agrees with the MC codes to the extent they agree with each other. In the present report, the 3DHZETRN code is extended to enable analysis in general combinatorial geometry. A more complex shielding structure with internal parts surrounding a tissue sphere is considered and compared against MC simulations. It is shown that even in the more complex geometry, 3DHZETRN agrees well with the MC codes and maintains a high degree of computational efficiency.

3DHZETRN: Neutron leakage in finite objects.

The 3DHZETRN formalism was recently developed as an extension to HZETRN with an emphasis on 3D corrections for neutrons and light ions. Comparisons to Monte Carlo (MC) simulations were used to verify the 3DHZETRN methodology in slab and spherical geometry, and it was shown that 3DHZETRN agrees with MC codes to the degree that various MC codes agree among themselves. One limitation of such comparisons is that all of the codes (3DHZETRN and three MC codes) utilize different nuclear models/databases; additionally, using a common nuclear model is impractical due to the complexity of the software. It is therefore difficult to ascertain if observed discrepancies are caused by transport code approximations or nuclear model differences. In particular, an important simplification in the 3DHZETRN formalism assumes that neutron production cross sections can be represented as the sum of forward and isotropic components, where the forward component is subsequently solved within the straight-ahead approximation. In the present report, previous transport model results in specific geometries are combined with additional results in related geometries to study neutron leakage using the Webber 1956 solar particle event as a source boundary condition. A ratio is defined to quantify the fractional neutron leakage at a point in a finite object relative to a semi-infinite slab geometry. Using the leakage ratio removes some of the dependence on the magnitude of the neutron production and clarifies the effects of angular scattering and absorption with regard to differences between the models. Discussion is given regarding observed differences between the MC codes and conclusions drawn about the need for further code development. Although the current version of 3DHZETRN is reasonably accurate compared to MC simulations, this study shows that improved leakage estimates can be obtained by replacing the isotropic/straight-ahead approximation with more detailed descriptions.

3DHZETRN: Shielded ICRU spherical phantom.

A computationally efficient 3DHZETRN code capable of simulating High (H) Charge (Z) and Energy (HZE) and light ions (including neutrons) under space-like boundary conditions with enhanced neutron and light ion propagation was recently developed for a simple homogeneous shield object. Monte Carlo benchmarks were used to verify the methodology in slab and spherical geometry, and the 3D corrections were shown to provide significant improvement over the straight-ahead approximation in some cases. In the present report, the new algorithms with well-defined convergence criteria are extended to inhomogeneous media within a shielded tissue slab and a shielded tissue sphere and tested against Monte Carlo simulation to verify the solution methods. The 3D corrections are again found to more accurately describe the neutron and light ion fluence spectra as compared to the straight-ahead approximation. These computationally efficient methods provide a basis for software capable of space shield analysis and optimization.

A comparative study of space radiation organ doses and associated cancer risks using PHITS and HZETRN.

NASA currently uses one-dimensional deterministic transport to generate values of the organ dose equivalent needed to calculate stochastic radiation risk following crew space exposures. In this study, organ absorbed doses and dose equivalents are calculated for 50th percentile male and female astronaut phantoms using both the NASA High Charge and Energy Transport Code to perform one-dimensional deterministic transport and the Particle and Heavy Ion Transport Code System to perform three-dimensional Monte Carlo transport. Two measures of radiation risk, effective dose and risk of exposure-induced death (REID) are calculated using the organ dose equivalents resulting from the two methods of radiation transport. For the space radiation environments and simplified shielding configurations considered, small differences (<8%) in the effective dose and REID are found. However, for the galactic cosmic ray (GCR) boundary condition, compensating errors are observed, indicating that comparisons between the integral measurements of complex radiation environments and code calculations can be misleading. Code-to-code benchmarks allow for the comparison of differential quantities, such as secondary particle differential fluence, to provide insight into differences observed in integral quantities for particular components of the GCR spectrum.

Dosimetric impacts of microgravity: an analysis of 5th, 50th and 95th percentile male and female astronauts.

Computational phantoms serve an important role in organ dosimetry and risk assessment performed at the National Aeronautics and Space Administration (NASA). A previous study investigated the impact on organ dose equivalents and effective doses from the use of the University of Florida hybrid adult male (UFHADM) and adult female (UFHADF) phantoms at differing height and weight percentiles versus those given by the two existing NASA phantoms, the computerized anatomical man (CAM) and female (CAF) (Bahadori et al 2011 Phys. Med. Biol. 56 1671-94). In the present study, the UFHADM and UFHADF phantoms of different body sizes were further altered to incorporate the effects of microgravity. Body self-shielding distributions are generated using the voxel-based ray tracer (VoBRaT), and the results are combined with depth dose data from the NASA codes BRYNTRN and HZETRN to yield organ dose equivalents and their rates for a variety of space radiation environments. It is found that while organ dose equivalents are indeed altered by the physiological effects ofmicrogravity, the magnitude of the change in overall risk (indicated by the effective dose) is minimal for the spectra and simplified shielding configurations considered. The results also indicate, however, that UFHADMand UFHADF could be useful in designing dose reduction strategies through optimized positioning of an astronaut during encounters with solar particle events.

Response functions for computing absorbed dose to skeletal tissues from neutron irradiation.

Spongiosa in the adult human skeleton consists of three tissues-active marrow (AM), inactive marrow (IM) and trabecularized mineral bone (TB). AM is considered to be the target tissue for assessment of both long-term leukemia risk and acute marrow toxicity following radiation exposure. The total shallow marrow (TM(50)), defined as all tissues lying within the first 50 µm of the bone surfaces, is considered to be the radiation target tissue of relevance for radiogenic bone cancer induction. For irradiation by sources external to the body, kerma to homogeneous spongiosa has been used as a surrogate for absorbed dose to both of these tissues, as direct dose calculations are not possible using computational phantoms with homogenized spongiosa. Recent micro-CT imaging of a 40 year old male cadaver has allowed for the accurate modeling of the fine microscopic structure of spongiosa in many regions of the adult skeleton (Hough et al 2011 Phys. Med. Biol. 56 2309-46). This microstructure, along with associated masses and tissue compositions, was used to compute specific absorbed fraction (SAF) values for protons originating in axial and appendicular bone sites (Jokisch et al 2011 Phys. Med. Biol. 56 6857-72). These proton SAFs, bone masses, tissue compositions and proton production cross sections, were subsequently used to construct neutron dose-response functions (DRFs) for both AM and TM(50) targets in each bone of the reference adult male. Kerma conditions were assumed for other resultant charged particles. For comparison, AM, TM(50) and spongiosa kerma coefficients were also calculated. At low incident neutron energies, AM kerma coefficients for neutrons correlate well with values of the AM DRF, while total marrow (TM) kerma coefficients correlate well with values of the TM(50) DRF. At high incident neutron energies, all kerma coefficients and DRFs tend to converge as charged-particle equilibrium is established across the bone site. In the range of 10 eV to 100 MeV, substantial differences are observed among the kerma coefficients and DRF. As a result, it is recommended that the AM kerma coefficient be used to estimate the AM DRF, and that the TM kerma coefficient be used to estimate the TM(50) DRF below 10 eV. Between 10 eV and 100 MeV, the appropriate DRF should be used as presented in this study. Above 100 MeV, spongiosa kerma coefficients apply well for estimating skeletal tissue doses. DRF values for each bone site as a function of energy are provided in an electronic annex to this article available at http://stacks.iop.org/0031-9155/56/6873/mmedia.

Response functions for computing absorbed dose to skeletal tissues from photon irradiation--an update.

A comprehensive set of photon fluence-to-dose response functions (DRFs) is presented for two radiosensitive skeletal tissues-active and total shallow marrow-within 15 and 32 bone sites, respectively, of the ICRP reference adult male. The functions were developed using fractional skeletal masses and associated electron-absorbed fractions as reported for the UF hybrid adult male phantom, which in turn is based upon micro-CT images of trabecular spongiosa taken from a 40 year male cadaver. The new DRFs expand upon both the original set of seven functions produced in 1985, and a 2007 update calculated under the assumption of secondary electron escape from spongiosa. In this study, it is assumed that photon irradiation of the skeleton will yield charged particle equilibrium across all spongiosa regions at energies exceeding 200 keV. Kerma coefficients for active marrow, inactive marrow, trabecular bone and spongiosa at higher energies are calculated using the DRF algorithm setting the electron-absorbed fraction for self-irradiation to unity. By comparing kerma coefficients and DRF functions, dose enhancement factors and mass energy-absorption coefficient (MEAC) ratios for active marrow to spongiosa were derived. These MEAC ratios compared well with those provided by the NIST Physical Reference Data Library (mean difference of 0.8%), and the dose enhancement factors for active marrow compared favorably with values calculated in the well-known study published by King and Spiers (1985 Br. J. Radiol. 58 345-56) (mean absolute difference of 1.9 percentage points). Additionally, dose enhancement factors for active marrow were shown to correlate well with the shallow marrow volume fraction (R(2) = 0.91). Dose enhancement factors for the total shallow marrow were also calculated for 32 bone sites representing the first such derivation for this target tissue.

The effect of anatomical modeling on space radiation dose estimates: a comparison of doses for NASA phantoms and the 5th, 50th, and 95th percentile male and female astronauts.

The National Aeronautics and Space Administration (NASA) performs organ dosimetry and risk assessment for astronauts using model-normalized measurements of the radiation fields encountered in space. To determine the radiation fields in an organ or tissue of interest, particle transport calculations are performed using self-shielding distributions generated with the computer program CAMERA to represent the human body. CAMERA mathematically traces linear rays (or path lengths) through the computerized anatomical man (CAM) phantom, a computational stylized model developed in the early 1970s with organ and body profiles modeled using solid shapes and scaled to represent the body morphometry of the 1950 50th percentile (PCTL) Air Force male. With the increasing use of voxel phantoms in medical and health physics, a conversion from a mathematical-based to a voxel-based ray-tracing algorithm is warranted. In this study, the voxel-based ray tracer (VoBRaT) is introduced to ray trace voxel phantoms using a modified version of the algorithm first proposed by Siddon (1985 Med. Phys. 12 252-5). After validation, VoBRAT is used to evaluate variations in body self-shielding distributions for NASA phantoms and six University of Florida (UF) hybrid phantoms, scaled to represent the 5th, 50th, and 95th PCTL male and female astronaut body morphometries, which have changed considerably since the inception of CAM. These body self-shielding distributions are used to generate organ dose equivalents and effective doses for five commonly evaluated space radiation environments. It is found that dosimetric differences among the phantoms are greatest for soft radiation spectra and light vehicular shielding.